Mean-Field Theory of Feshbach-Resonant Interactions in 85 Rb Condensates 
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Recent Feshbach-resonance experiments with 85 Rb Bose-Einstein condensates have led to a host of 
unexplained results: dramatic losses of condensate atoms for an across-resonance sweep of the mag- 
netic field, a collapsing condensate with a burst of atoms emanating from the remnant condensate, 
increased losses for decreasing interaction times — until very short times are reached, and coherent 
oscillations between remnant and burst atoms. Using a simple yet realistic mean-field model, wc find 
that rogue dissociation, molecular dissociation to noncondensate atom pairs, is strongly implicated 
as the physical mechanism responsible for these observations. 
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Introduction- The process known as the Feshbach res- 
onance occurs when two ultracold atoms collide in the 
presence of a magnetic field, whereby a spin flip of one 
atom can induce the pair to jump from the two-atom 
continuum to a quasibound molecular state. If the ini- 
tial atoms are Bose condensed, the so-formed molecules 
will also comprise a Bose-Einstein condensate (BEC) [0. 
Since the Feshbach resonance is mathematically identi- 
cal to photoassociation the process that occurs 
when two ultracold atoms form a molecule by absorb- 
ing a photon, insight gathered in either case is applicable 
to the other. In particular, the most recent results from 
photoassociation predict that rogue dissociation, molec- 
ular dissociation to noncondensate atom pairs, imposes a 
maximum achievable rate on atom-molecule conversion, 
as well as the possibility of coherent Rabi oscillations be- 
tween the BEC and the dissociated atom pairs || . 

As predicted ||, the Feshbach resonance was used to 
enable Bose condensation in 85 Rb by tuning the natively 
negative atom-atom scattering length into the positive 
regime Rj. Part of this experiment involved a sweep of 
the magnetic field across the Feshbach resonance, result- 
ing in heavy condensate losses (~ 80% for the slowest 
sweep rates). Additional experiments led to the obser- 
vation of collapsing condensates, an event characterized 
by bursts of atoms emanating from a remnant BEC, and 
coined "bosenova" for the analogy with a supernova ex- 
plosion Q . More recently, experiments with pulsed mag- 
netic fields that come close to, but do not cross, the 
Feshbach-resonance have revealed a striking increase in 
condensate loss for a decrease in the interaction time- un- 
til reaching very short times |J . Finally, double-pulse re- 
sults indicate large amplitude remnant-burst oscillations 
(~ 25%), with the missing atoms prompting speculation 
on the formation of a molecular condensate (lO) . 

However, recent work on Feshbach-stimulated photo- 
production of stable molecular condensates indicates that 
rogue dissociation dominates atom-molecule conversion 
for the above 85 Rb Feshbach resonance, meaning that it is 
practically useless for producing any significant fraction 



of molecular BEC, stable or otherwise JT1J. Additionally, 
assertions of a breakdown of mean-field theory in the face 
of large resonance-induced scattering lengths |?],|| are in- 
triguing, especially given that theory actually faults the 
effective all-atom description Jl^]. The purpose of this 
Letter is to test a mean-field model of coherent atom- 
molecule conversion against the salient features of the 
JILA experiments, with the intent of developing an un- 
derstanding of the basic physics involved. 

Mean-field theory- The mathematical equivalence of 
the Feshbach-resonance [|J and photoassociation ||,|| lies 
in both processes being described, in terms of second 
quantization, as destroying two atoms and creating a 
molecule. We therefore model a quantum degenerate gas 
of atoms coupled via a Feshbach resonance to a conden- 
sate of quasibound molecules based on Refs. The 
initial atoms are denoted by the boson field 4>(r,t), and 
the quasibound molecules by the field ip(r, t). The Hamil- 
tonian density for this system is 
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where m — 2/i is the mass of an atom, Mo is the energy 
difference between a molecule and two atoms, p is an 
invariant density equal to the sum of atom density and 
twice the molecular density, and T(e) is the dissociation 
rate for a molecule with the energy fie above the threshold 
of the Feshbach resonance. 

Switching to momentum space, only zero-momentum 
atomic and molecular condensate modes are retained, 
represented by the respective c-number amplitudes a and 
(3. We also take into account correlated pairs of noncon- 
densate atoms using a complex amplitude C(e), which 
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represent pairs of noncondensate atoms in the manner of 
the Heisenberg picture expectation value (a p a_ p ), with 
Tie being the relative energy of the atoms. The nor- 
malization of our mean fields is such that \a\ 2 + \(3\ 2 + 
J de\C(e)\ 2 = 1. We work from the Heisenberg equa- 
tion of motion of the boson operators under the simplify- 
ing assumption that the noncondensate atoms pairs are 
only allowed to couple back to the molecular conden- 
sate, ignoring the possibility that noncondensate atoms 
associate to make noncondensate molecules. This neglect 
is justified to the extent that Bose enhancement favors 
transitions back to the molecular condensate. The final 
mean-field equations are H 
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The analog of the Rabi frequency for the rogue modes 
£ is inferred using Fermi Golden rule, which gives the 
dissociation rate for a positive-energy molecule as 
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Next the problem is reformulated in terms of two 
key parameters with the dimension of frequency. The 
density-dependent frequency 
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has been identified before, along with the operational 



significance that, once 51 



rogue dissociation is ex- 



pected to be a dominant factor in the dynamics |j] |). 
Here it is convenient to define another primary parame- 
ter with the dimension of frequency. Considering on-shell 
dissociation of molecules to atoms with the relative en- 
ergy e, the Wigner threshold law delivers a dissociation 
rate T(e) such that T(e)/y/e converges to a finite limit for 
e — > 0; hence, we define 
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which indeed has the dimension of frequency. Compar- 
ison of Eqs. (g), (^), and (^3j) give the parameters in 
the mean-field equations as 
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Renormalization- When the coupling to the continuum 
of noncondensate atom pairs is included, the continuum 
shifts the molecular state [Q. As explained before ||, 
and will be discussed in more detail below, we have taken 
the dominant state pushing into account in our calcu- 
lations. However, there is a relevant residual effect to 
consider. 



To begin with, consider the equation of motion for the 
molecular amplitude, including the coupling to noncon- 
densate atom pairs but not to the atomic condensate [set 
51 = in Eqs. (^)j. Again, the particular energy depen- 
dence of the coupling comes from the Wigner threshold 
law, which is here assumed to be valid for all relevant 
energies; but, as is necessary in the numerical calcula- 
tions anyway, we cut off the coupling between molecules 
and atom pairs at some frequency eu- The question of 
renormalization is the question of the dependence of the 
results on the cutoff e^f. 

Equations @ (with 51 = 0) are easy to solve using, say, 
the Fourier transformation. With the initial condition 
that Pit = 0) = 1, for positive times the solution has the 
Fourier transform 



cu + ir) — Sq — S(w) 
where the self-energy is 
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and r\ = + . Now, a real pole of /3(w) corresponds to 
a true stationary state of the Hamiltonian. It turns 
out that for suitable detunings there is a real pole 
lu < 0, which obviously corresponds to the coupling- 
renormalized energy of the molecules. 

Assuming that the value of to is negative, the integral 
@ is carried out easily and the equation giving the pole 
becomes 
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The term involving the detuning is the main contribution 
to the renormalization. As we have done before in our 
numerical calculations ||, we choose the bare detuning 
So so that, for the given energy cutoff ejw, the renormal- 
ized detuning So — s/cmS/tt attains the desired value. 
Hereafter we use the symbol S for the renormalized de- 
tuning. This is the parameter that is varied by changing 
the laser frequency in photoassociation, or the magnetic 
field in the Feshbach resonance. We carry out the rest 
of the renormalization by setting cm — > oo. We find the 
equation giving the characteristic frequency correspond- 
ing to the molecules as lu — 5 — V — = . The proper 
solution is 
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(ii) 



which is valid for all negative detunings and gives the 
characteristic frequency of the molecule. Nevertheless, 
complete diagonalization of the problem (with 51 = 0) 
shows that the mode evolving at this frequency is not 
the original or "bare" molecules, but a coherent super- 
position of molecules and noncondensate atoms pairs. 
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Numerical procedures- The mean-field equations (|3j) 
are integrated using the norm-preserving predictor- 
corrector algorithm described in Ref . || . Magnetic fields 
are converted to detunings according to hS = A[j,(Bq — 
B), where the position of the Feshbach resonance is 
So = 154.9 G, and where the difference in magnetic mo- 
ments between bound molecules and free atom pairs, 
A M w 2 fiB (where /is is the Bohr magneton), is bor- 
rowed from 87 Rb |Q. Compared to the ensuing detun- 
ings 5, the interactions energies between the atoms due to 
the background scattering length a — 23.8 nm are imma- 
terial. We therefore ignore atom-atom interactions unre- 
lated to the Feshbach resonance, as well as the (unknown) 
atom-molecule and molecule-molecule interactions. 

We have used a number of different methods to es- 
timate the parameter S, all giving similar values. The 
present argument goes as follows. One of the experi- 
ments in Ref. jlO| gives the characteristic frequency (pre- 
sumably) of molecules, u) = —2.07 x 10 5 s _1 , at the mag- 
netic field B = 159.69 G corresponding to the detuning 
5 = -8.42 x 10 7 s _1 . Solving from Eq. (O), we have 
5 = 5.29 x 10 9 s-\ and thus £ = 381 s" 1 ^ It should 
be noted that, while our rogue-dissociation coupling can 
be shown to give the correct atom-atom scattering the- 
ory close to the resonance (5 — ► 0), Eq. ( p"lj ) will not 
correctly reproduce the molecular energy on the side of 
large magnetic fields in Fig. 5 of Ref. JTo|] . Passable agree- 
ment could be reached by treating the energy cutoff tu 
as a variable finite parameter; unfortunately, the abrupt 
cutoff at 6m introduces physical and numerical artifacts, 
which need to be cleaned up via a judicious energy de- 
pendence of the coupling. For example, the entire exper- 
imental range of magnetic fields can be fit by multiplying 
the original coupling by an exponential [exp(— e/em)] or 
rational [{^/(e + fM) 2 ] cutoff. As our aims are strictly 
qualitative, we discuss these details elsewhere. 

Results- We begin with the Cornish et al. experi- 
ments implementing a sweep of the magnetic field across 
the Feshbach resonance, which is of course a version of 
the age-old Landau-Zener problem f^jUJIs}]. Although a 
sweep of the detuning S from above to below threshold 
at a rate slow compared to the condensate coupling fl 
will move the system adiabatically from all atoms to 
all molecules, rogue dissociation will overtake coherent 
atom-molecule conversion when f2 ~ ui p Never- 
mind that the JILA experiments sweep from below to 
above threshold, for a density p — 1 x 10 12 cm -3 the con- 
densate coupling is SI = 1.93 x 10 5 s~ 1 « 250 ui p , and 
so rogue dissociation should seriously dominate. This is 
indeed the case (see Fig. [I]). Apparently, coherent conver- 
sion occurs not between atomic and molecular BEC, but 
between atomic BEC and dissociated atom pairs. Hold- 
ing this thought, we conclude that mean-field theory in- 
dicates rogue dissociation as a primary sink of atoms in 
the Ref. M sweeps across the Feshbach resonance. 

Next we consider the experiments of Claussen et al. , 
for which nontrivial electromagnetic coil technology was 



developed to create trapezoidal magnetic field pulses that 
bring the system near- but not across- resonance, hold 
for a given amount of time, and return to the original 
field value. Neglecting the burst, these remnant-focused 
experiments revealed a contradiction with the conven- 
tional understanding of condensate loss: rather than a 
loss that increased monotonically with increasing inter- 
action time, the results indicated a loss that increased 
with decreasing interaction time, until very short times 
were reached. The present mean-field approach works 
similarly, as shown in Fig. ||. Our interpretation is that 
adiabaticity is again at play. At very short pulse du- 
rations, increasing interaction time leads to increasing 
condensate loss, as expected. In contrast, as the time de- 
pendence of the pulse gets slower, the system eventually 
follows the pulse adiabatically, and returns close to the 
initial condensate state when the pulse has passed. 

Finally, we turn to the experiments performed by Don- 
ley et al. [ho|, in which two trapezoidal pulses were ap- 
plied to a ™Rb condensate, and the fraction of remnant 
and burst atoms measured for a variable between-pulse 
time and magnetic-field amplitude. These experiments 
revealed coherent remnant-burst oscillations with ampli- 
tudes of up to ~ 25%. As it happens, we have recently 
predicted coherent oscillations between atoms and dis- 
sociated atom pairs in a rogue-dominated system, al- 
though we harbored doubts regarding any practical re- 
alization ||. Casting these doubts aside, we consider 
a time dependent detuning (magnetic field) in essence 
lifted from Fig. 2 of Ref. Q [Fig. ||(a)], and determine 
the fraction of remnant condensate atoms, nonconden- 
sate atoms, and molecules at the end of the pulse se- 
quence as a function of the holding time between the two 
pulses [Fig. ||(b)]. Oscillations are seen with the ampli- 
tude of about 15% between condensate and nonconden- 
sate atoms at the frequency of the molecular state corre- 
sponding to the magnetic field during the holding period. 
The molecular fraction appears too small to account for 
the amplitude of the oscillations. In fact, what we termed 
molecular frequency is the characteristic frequency of a 
coherent superposition of molecules and noncondensate 
atom pair. Here the oscillations, directly comparable to 
Fig. 4(a) in Ref. Jl(J , are Ramsey fringes |ll| in the evo- 
lution between an atomic condensate and a molecular 
condensate dressed with noncondensate atom pairs. 

Although our rogue-dissociation ideas provide a neat 
qualitative explanation for the three experiments we have 
discussed jTfl], in all fairness it must be noted that we 
have fallen short of a full quantitative agreement. It ap- 
pears that our model is missing a so far unidentified ad- 
ditional loss mechanism for the condensate atoms. 

Conclusions- We have demonstrated that a minimal 
mean-field model is sufficient to qualitatively explain a 
number of puzzling results in Feshbach-resonant systems. 
On the whole, collapsing-condensate physics is therefore 
understood as a matter of rogue dissociation, which leads 
to strong losses in the threshold neighborhood, decreased 
remnant fraction for decreasing interaction time — until 
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very short times are reached, and coherent remnant-burst 
oscillations. Furthermore, although atom-molecule co- 
herence has no doubt been achieved Jul, the amplitude 
of the remnant-burst oscillations need not be indicative 
of the number of condensate molecules present. Ironi- 
cally, the strength of the Feshbach resonance has led to a 
regime dominated by rogue dissociation, which tends to 
counteract the production of a molecular condensate. 
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FIG. 1. (a) Experimental [ij (□) and theoretical (o) atom 
loss incurred in sweeping a D Rb BEC across the Feshbach 
resonance, where the magnetic field is swept in a linear fash- 
ion from Bi = 162 G to Bf — 132 G. In each numerical run, 
the fraction of molecular condensate is ~ 10 -6 . (b) Results 
for S _1 = 100 us/G are typical, and suggest that the system 
undergoes collective adiabatic following from BEC to dissoci- 
ated atom pairs. 
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FIG. 2. Theory of a magnetic field pulse applied to a 85 Rb 
condensate for p = 1.9 x 10 13 cm" 3 and Q. = 8.42 x 10 5 s" 1 . 
(a) Remnant fraction versus detuning (magnetic field) rise 
time, (b-d) Results for a pulse with 150 fjs rise time indicate 
adiabatic passage of BEC atoms to and from dissociated atom 
pairs. The minimum in panel (a), similar to Ref. 0], occurs 
at the onset of adiabaticity. 
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FIG. 3. Simulation of the Ref. [jlo| experiments for a den- 



sity p = 5.4 x 10 13 cm" 3 and Q = 1.42 x 10 6 s"\ (a) Time de- 
pendence of the detuning, and (b) the fraction of atoms in the 
remnant condensate (solid line) , in noncondensate atoms pairs 
(dashed line) and in the molecular condensate (short-dashed 
line) after the pulse sequence as a function of the hold time 
th between the two trapezoidal pulses. The frequency of the 
oscillations is compatible with the prediction from Eq. (pd[), 
identifying these oscillations as Ramsey fringes in the transi- 
tion between the atomic condensate and a molecular conden- 
sate dressed by noncondensate atom pairs. 
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